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ABSTRACT 



We present results from the Radio Interferometric Planet (RIPL) search for compan- 
ions to the nearby star GJ 896A. We present 11 observations over 4.9 years. Fitting as- 
trometric parameters to the data reveals a residual with peak-to-peak amplitude of ~ 3 
mas in right ascension. This residual is well-fit by an acceleration term of 0.458 ± 0.032 
mas/y 2 . The parallax is fit to an accuracy of 0.2 mas and the proper motion terms are 
*~^ fit to accuracies of 0.01 mas/y. After fitting astrometric and acceleration terms resid- 

uals are 0.26 mas in each coordinate, demonstrating that stellar jitter does not limit 
the ability to carry out radio astrometric planet detection and characterization. The 
acceleration term originates in part from the companion GJ 896B but the amplitude 
of the acceleration in declination is not accurately predicted by the orbital model. The 
acceleration sets a mass upper limit of 0.15 Mj at a semi-major axis of 2 AU for a 
planetary companion to GJ 896A. For semi-major axes between 0.3 and 2 AU upper 
limits are determined by the maximum angular separation; the upper limits scale from 
the minimum value in proportion to the inverse of the radius. Upper limits at larger 
r**'» ■ radii are set by the acceleration and scale as the radius squared. An improved solution 

for the stellar binary system could improve the exoplanet mass sensitivity by an order 
of magnitude. 

Subject headings: astrometry, stars: activity, stars: early-type, stars: planetary systems, 
H ! radio continuum: stars 
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Introduction 



The study of extrasolar planets provides an important link between the study of star forma- 
tion, the study of our own solar system, and the search f or extraterrestrial li fe. Radial velocity 
searches 



Butler et al. 



searches (Gaudi et al 



2006), transit searches with Kepler (JBorucki et al.ll201ll ). and microlensing 



200a) have surveyed significant components of the exoplanet mass-separation 



parameter space exposing a wide variety of phenomenology. 

Despite their abundance, low mass stars have been less well-characterized for a number of 
reasons. In particular, low mass stars are often too faint or too strongly variable for high precision 
radial velocity and transit searches. Microlensing searches have ide ntified companions to low mass 



stars (Bennett et al. 



2008 



Beaulieu et al. 



2006 



Gould et al.l 120061 ) . which suggest that such com- 



panions may be very common. On the oth er hand, radial ve l ocity searches sugges t very low planet 



occurrence rates for close-in companions (JEndl et al. 



2006 



Johnson et al.l 120071 ) . The frequency 



and distribution of companions to low mass stars can be an important diagnostic of m ethods for 
planet formation (JLaughlin et alJbood : Eda & Li J 12009 : IbossI 12006 : kennedv et al.ll2007l ). 



Astrometry ca n be a powerful tool for detection and characterization of nearby, low-mass 
planetary systems. iMuterspaugh et al.l (|2010l ) demonstrate that optical astrometric measurements 
with the Palomar Testbed Interferometer are capable of achieving sub-milliarcsecond accuracy in the 
separation between close binaries. They summarize the known properties of sub-stellar companions 
to low mass stars in binary systems and conclude that giant pl anet companions in such systems are 
common. The Carnegie Astrometric Planet Search Program (JBoss et al.ll2009l ) has also obtained 
sub-milliarcsecond astrometry on the low mass star NLTT 48256 over 2 years of observations; this 
program will ultimately survey 100 nearby M, L, and T dwarfs. 

The Very Long Baseline Array (VLBA) can routinely achieve an astrometric accuracy of ~ 100 
//as in a single epoch at frequencies ~ 10 GHz and ~ 20 //as at high er frequencies. The VLBA is 
capable of accuracies as high as 8 //as under favorable circumstances (JFomalont fe Kopeikinll2003l ). 
To use this technique the target source must have a sufficiently high brightness temperature to be 
detected by a high resolution radio interferometer. For the VLBA, brightness temperatures must 
be T;, ~ 10 7 K, requiring nonthermal emission. Thus, the systems which can be studied are limited 
to the most active stars. VLBA a strometry has been important for a wide r ange of stellar mea- 



surements includin g parallax (e.g.. ISandstrom et al 



2007 



Menten et al. 



2011 



2007; lLoinard et all 120071 : 



Brunthaler et al. 



2011 



Torres et al.l 12009 ), galactic structure (|Xu et al.l 120061 . _ 

chara cterization (Guira do et al.l 1200m. an d the radio/optical reference- frame link ( Lestrade et al 

1999J). Recently, Forbrich & Berger (12009!) detected a brow n dwarf with the VLBA and explored 



, binary orbit 



the possibilities of companion detection. iBower et al.l (J2007I ) reviewed the prospects for future radio 
astrometric surveys with respect to other methods. 

RIPL, the Radio Interferometric Planet Search, is a 4- year program using the VLBA and 
the 100m Green Bank Telescope (GBT) to search for the astrometic signatures of massive planets 
around nearby {D < 10 pc), low-mass (M dwarf) stars. We are observing a sample of 30 stars a total 



of 12 times over a 4-year period with the expectation of achieving ~ 0.1 milliarcsecond astrometric 
accuracy per epoch, and ultimately having the sensitivity to detect Jupiter-mass planets at 1 AU. 
This astrometric search is an important complement to radial velocity, transit, microlensing, and 
optical astrometric searches for planets on account of greater sensitivity to 1) low mass stars, 2) 
long-period planets, 3) active stars, and, 4) planets that may be imaged with extreme adaptive 
optics. Many RIPL stars are in binary systems, which may have higher planetary fractions. 



In Paper I (JBower et al.ll2009l ). we used Very Large Array (VLA) observations to identify the 
sample of stars appropriate for RIPL and we used the VLBA to characterize the detectability and 
short-term (~ 10 days) stability of the astrometric images of these stars. The short-term stability 
addresses a significant concern of radio astrometry that stellar activity will lead to large astrometric 
errors. The absence of a significant difference in the historical optical proper motions and the radio 
proper motion from this RIPL precursor survey permitted us to exclude brown dwarf companions 
to the four stars observed. RIPL observing began in October 2007 and is anticipated to end in late 
2011. 

In this paper, we present initial results from RIPL for the star GJ 896A. This star was also 
observed as part of the precursor survey in paper I. Together these observations demonstrate 
astrometric stability over a 4.9 year span. GJ 896A (EQ Peg A) is a variable star in a multiple 
system with GJ 896 B, C, and D. GJ 896B i s its closest companion; the orbit has been characterized, 
but it is not well constrained ()Heintzlll984l ). GJ 896 A and B have spectral types M3.5 and M4.5. 
These two stars have bee n detected previously in the radio and seen t o have point-like s tructure and 



stron g radio variability (jPallavicini et al.l Il985l : iJackson et al.1 Il989i : iBenz et al.1 Il995l : iGagne et al. 



1998). In §2, we describe our observations, analysis, and results. In §3, we present our astrometric 
analysis of the data. In §4, we examine the constraints that we can place on companions to GJ 
896A. In S5, we summarize our results. 



2. Observations, Data Analysis, and Results 

Observations of GJ 896A were obtained as part of the RIPL survey. RIPL observations spanned 
October 2007 through June 2011, and are ongoing. In each epoch, we observed two stars over a span 
of 8 hours. Fast switching with a phase reference calibrator and secondary (and sometimes tertiary) 
astrometric calibrators reduced total observing time on each star to approximately two hours. The 
total span from start to end of observation for a particular star within an epoch varied, depending 
on the relative coordinates of the stars and the allocated time slot. In some instances, the stars 
were observed sequentially in four hour blocks. In other instances, the stars were interleaved and 
the span was over the full eight hours. 

Observations typically included all ten antennas of the Very Long Baseline Array as well as 
the 100 meter Green Bank telescope. Data were recorded at a rate of 512 Mbps in dual circular 
polarization mode with observing bands at 8.4 GHz using the standard VLBA mode v4cm-512- 



8-2. More precisely, this mode corresponds to four dual-polarization frequency bands centered 
at 8384.49, 8400.49, 8416.49, and 8432.49 MHz, each sampled at 32 Msamples/second with 2 bit 
sampling. Occasional observations of fringe calibrators were interspersed throughout each epoch. 
The switching cycle included repeated observations of 1 minute on the phase reference calibrator 
followed by 3 minutes on the star. Every 30 minutes short duration observations of the secondary 
and tertiary calibrators phase referenced to the primary calibrator were obtained. 



Data were reduced with an AIPS (|Greisenll2003l ) pipeline that performed standard calibrations 



of amplitude and phase. Corrections for the atmospheric total electron content (TEC) and post- 
correlation updates to the Earth orientation parameters (EOP) were also included. All sources 
were imaged and deconcolved using the CLEAN algorithm. The typical synthesized beam for GJ 
896A was 3.0 x 1.1 mas, oriented North-South. We fit Gaussian models to all sources using the 
AIPS task JMFIT. 

In the case of GJ 896A, the primary phase reference calibrator is J2328+1929, which has an 
ICRF2 position that is accurate to 1.84 mas x 1.58 masllj. The assumed position for J2328+1929 is 
(23 h 28 m 24.874773 s , 19° 29'58.03010"); these updated values are slightly offset from the assumed 
position of paper I by (0.254, 0.059) mas in right ascension and declination, respectively. All 
positions for other sources are relative to this position. The secondary calibrator is J2334+2010. 
No tertiary calibrator was used. We apply a correction for tropospheric phase gradients on the 
astrometric position of the star by removing variations in the position of the secondary calibrator. 
This correction is applied on an epoch-by-epoch basis. A more sophisticated approach using the 
AIPS task ATMCA can make corrections on shorter timescales. The magnitude of corrections to GJ 
896 A are ~ 0.1 mas. This approach may overestimate trophospheric errors slightly given the greater 
distance and higher declination of J2334+2010 from the primary calibrator. The close spacing and 
nearly linear arrangement of the star and calibrators, however, make this a good calibration set 
(Fig. [TJ). Typically, 3C 454.3 was observed as the fringe calibrator. As of June 2011, we have 8 
good epochs for GJ 896A. 

Both calibrators are dominated by a strong point source and show only weak evidence for 
non-point-like structure. Images of the calibrators are included in Paper I. The primary calibrator, 
J2328+1929, has a typical flux density of 60 mJy; the secondary calibrator, J2334+2010, has a 
typical flux density of 20 mJy. In Table [TJ we summarize flux densities and positions for the 
secondary calibrator as a function of epoch; we also include the mean value from the three epochs 
of paper I. The positions for J2334+2010 have an rms variation of (0.11, 0.09) mas and a mean 
offset (0.06, 0.22) mas from the astrometric position obtained in paper I. That mean offset is within 
the error bars of the original measurement. These measures for the secondary calibrator provide 
an estimate of the astrometric accuracy of the experiment. 

We summarize positions and flux densities for G J 896A in Table [2] as a function of epoch. We 



: http: //gemini .gsfc .nasa.gov/solutions/2010a/2010a.html 



also include measurements from Paper I. GJ 896A is strongly variable: flux densities vary from 
0.140 to 9.512 mJy during the course of the experiment. The images of GJ 896A are predominantly 
point-like, as found in Paper I. The smaller errors on the flux density and the positions are reflective 
of the increased sensitivity and (u, v) coverage of RIPL over the precursor experiment due to the 
doubling of the observing time, doubling of the bit rate, and addition of the GBT to the array. The 
statistical errors are significantly less than the systematic errors associated with the astrometry. 
For the analysis discussed below, we add in quadrature an error of 0.2 mas to each coordinate as 
an estimate of the systematic error. 



3. Astrometric Fitting 

In addition to the eight epochs from RIPL, we obtained three epochs in the RIPL precursor 
survey over three days in March 2006. We consider the data from all 11 epochs in our astrometric 
fitting here. The data span 4.9 years. 

GJ 896A was not originally part of th e Hipparcos catalog but it is included in the reanalysis 



of the Hipparcos data (Ivan Leeuwenl 120071 ) . We summarize the known astrometric parameters in 



Table El GJ 896A is also known to be in a double sys tem with GJ 896B; parameters ar e determined 



from ~ 40 years of optical observations (Tabled iHeinta Il984l ; iMason et al.ll200ll ). Since, the 
baseline for this orbit determination is ~ 0.1 times the estimated period, the orbital parameters 
are poorly determined. This is reflected in the binary orbit grade of 5 ("indeterminate") on a scale 
of 1 to 5. 

In Table El we show the nominal astrometric data available from optical sources as well as best- 
fit solutions based on RIPL data. We performed two kinds of fits: one with just the astrometric 
parameters, and the other with astrometric parameters and an acceleration term. Figures El El 
and [H show the astrometric data, the models, and the residuals for the cases of the initial model, 
the astrometric model, and the astrometric plus acceleration model, respectively. The astrometric 
plus acceleration model that we fit to our data (a(£j),e)(£j)) at epoch tj is 

a(ti) = a + fi a (ti-t ) + Tr a (a,5,ti) + a a (ti-tQ) 2 , (1) 

S(ti) = <5o + ^(ti-io) + ^(a,<5,t i ) + a (5 (i,-4) 2 , (2) 

where (ao,5o) is the equinox J2000 position at time to, the epoch 2000.0, ({J. a ,fJ-s) is the proper 
motion, (ir a ,irs) are the projections of the parallax it onto the two coordinates, (a Q , ay) is the local 
acceleration, and t$ and t$ are zero-points for the acceleration in each coordinate. The initial model 
clearly shows a substantial error in \x a and a ~ 1 mas error in ir. Fitting the astrometric parameters 
(ao, <5o, fJ-a, fJ-8, tt) substantially improves the quality of the fit but leaves a parabolic residual in 
a. An acceleration term reduces the rms in a from 1.2 to 0.26 mas. The residual in declination is 
improved from 0.33 to 0.26 mas through the addition of an acceleration term. The astrometric fit 
alone with 22 measurements and 5 parameters producing 17 degrees of freedom results in x 2 = 338. 



The acceleration term significantly reduces x 2 to 29.7 with 9 parameters producing 13 degrees of 
freedom. The reduced x 2 = 2-3 in this final fit. 

We estimate errors in individual parameters by mapping the x 2 surface as a function of each 
parameter. For the purely astrometric case, the x 2 curves are parabolic. We plot the x 2 curves for 
the astrometric and acceleration terms in Figure These curves are also parabolic and symmetric 
about their best-fit values. We estimate the la errors from these curves. The astrometric parameter 
estimates are significantly more accurate when the acceleration terms are included. We do not 
estimate the full covariance matrix of the parameter estimates due to the computational cost of 
exploring a 9-dimensional space. We did explore the covariance between the parallax and the right 
ascension acceleration term, two parameters that one might expect to be strongly linked. We find 
that the x 2 surface in this two-dimensional domain is well-characterized by a paraboloidal surface 
with axes that are misaligned from the principle axes by less than 30 degrees; that is, the parameters 
are essentially independent. 

The estimated values for the astrometric parameters are determined one to two orders of 
magnitude more accurately than the Hipparcos values. The values for ao, So, tt, and fis are in 
agreement with the Hipparcos values. The proper motion in right ascension, fi a , however, differs 
by 16.5 ± 1.4 mas/y. Over the time baseline between the Hipparcos epoch (1991.25) and the mean 
RIPL epoch (2008.75), this implies an acceleration of 0.94 ± 0.08 mas/y 2 . This has the same sign 
and a comparable magnitude to the acceleration term measured from the RIPL data alone. The 
agreement in fi$ implies an upper limit of 0.06 mas/y 2 over this interval. 

Motion of the star during the observing epoch can be an important effect, especially if the 
flux density is variable on short timescales. For GJ 896A, the median absolute value of motions 
from astrometric parameters are (0.10, 0.048) mas/hour in the two coordinates, respectively, and are 
always lower than (0.14, 0.065) mas/hour. This can lead to smearing out of the source and reduction 
of the peak flux density. If the source flux density is strongly variable during the observation, then 
the mean time of the epoch may not correspond to the time of the astrometric detection, leading 
to errors be as large as (0.6,0.3) mas. The small residuals we obtain from fitting indicate that this 
effect is not likely to be important in most epochs. Future reductions of this data will attempt to 
determine any light curve variability and its effect on the astrometry and/or include corrections to 
the visibilities for the assumed astrometric motion of the star. In high signal to noise ratio cases, 
we may be able to simultaneously fit a position and velocity within a given epoch. 



4. Constraints on Companions to GJ 896A 

The best fit accelerations for GJ 896A are (0.458, 0.087) mas/y 2 . We compare this to the 
known parameters of the binary system. We plot relative separation, velocity, and acceleration 
determined from the binary parameters (Figure [6]) . Note that the definition of the orbit gives the 
position of B relative to A. Thus, the sign of velocity and acceleration terms is inverted relative 



to our measurement. If the stars have equal mass, then the motion relative to the center of mass 
will be half of the value plotted. The mean calculated accelerations for A relative to B over our 
observing epochs is (0.3 ± 0.1, 3.1 ± 0.6) mas/y 2 . The calculated acceleration term in a agrees 
reasonably well with the measured value and is comparable to the mean acceleration over the 
entire orbit ~ 2irr/P 2 = 0.3 mas/y 2 . The calculated acceleration term ag does not agree at all with 
the measured value. The RIPL measurements occur near the predicted periastron, and, therefore, 
a maximum in the value of as- Our results clearly indicate a significant error in the model for this 
binary system and appear to indicate that the period may be longer than 359 years, or that the 
eccentricity or time of periastron have not been properly estimated. 

The acceleration a allows us to a place a limit on the stellar companion mass, Mb, using 

(^A = ( a \ (-Y (3) 

\Mj) VO-0381 AU/y 2 y VlAU/ ' VJ 

where r is the current separation projected onto the sky and Mj is the mass of Jupiter. If we 
translate the measured acceleration to physical units, we find a a = 2.8 x 10~ 3 AU/y 2 and ag = 
0.5 x 10 -3 AU/y 2 , or the total acceleration a = 2.9 x 10~ 3 AU/y 2 . Applied to the stellar companion 
at a projected distance of 5. 5" (the optically measured separation in 2007), we estimate a mass for 
GJ 896B of Mb ~ 0.4 Mm. This mass es timate is in rough agreement with expectations for a star 



of spectral type M4.5 (JKraus et al.1 12011 



We can also apply the acceleration as a limit to estimate planetary mass companions. If 
we take the acceleration a as an upper limit then we can place a limit on a planetary mass, 
M p ~ 0.08Mj (r/lAU) . This acceleration limit applies to companions with a period ~ 4.9 y, or 
at r ~ 1.7 AU for a stellar mass of 0.2 Mq. For closer orbits, we can set a constraint based on 
the maximal amplitude of angular displacement. We set this limit by computing the periodogram 
for the residuals in each coordinate after subtraction of the astrometric plus acceleration model 
(Figure [7j). The amplitude of the periodogram falls well below the 95%-confidence threshold; we 
find no evidence for a short-period companion. The maximal power in the periodogram corresponds 
to an upper limit of 0.4 mas in angular displacement. We are sensitive to orbital periods as short 
as the characteristic separation between epochs ~ 140 days; this limit translates to an inner radius 
cutoff of ~ 0.3 AU. In Figure [8j we combine the acceleration limits and the angular displacement 
limits. The maximal sensitivity is at 2 AU with an upper limit of 0.15 Mj. Thus, we are sensitive 
to a planet with approximately 3 Neptune masses at large radii from GJ 896A. Note that this 
limit on the mass applies strictly to face-on systems; orbits at all inclinations will have a maximum 
acceleration equal to the face-on case but will have smaller instantaneous values depending on the 
phase of the orbit. 

A caveat to this calculation is that the acceleration we measure may not be a upper limit 
to the effect of a planetary companion if the companion and GJ 896B produce accelerations that 
almost cancel each other. The agreement in a a between measured and calculated values suggests 
that this is not likely; the disagreement in as could be due to a planetary companion acceleration 
that wholly or partially cancels the binary acceleration. A Jupiter mass companion at 2 AU would 



produce the magnitude of acceleration necessary for the cancellation. This cancellation, however, 
requires a coincidence that is unlikely to persist over a period of observation of this duration. The 
agreement in fj,g between the Hipparcos and RIPL measurements indicates that such a cancellation 
must persist over a period of ~ 16 y, which is improbable. Improved optical astrometry of the 
binary system can determine the contribution of GJ 896B to the motion of GJ 896A. 

The rms residuals indicate that the centroid of the stellar emission is confined to a physical 
region ~ 0.3R(7). T his is comparable to or less than the photospheric radius of a late- type star (e.g., 



Kraus et al.ll201ll ). These residuals occur over a range of source flux densities that vary by a factor 
of ~ 60. We conclude that the magnetic activity associated with GJ 896A is symmetric about the 
stellar surface and/or confined primarily to the disk of the star. 



5. Conclusions 

We have presented here an analysis of RIPL observations of the star GJ 896A. We measure 
the parallax to an accuracy of 0.2 mas and the proper motion to an accuracy of 0.01 mas/y. We 
measure an acceleration term of 0.466 ±0.045 mas/y 2 . This term is significantly less than expected 
based on the model for the orbit of the GJ 896AB system, suggesting an error in the optically- 
determined orbital parameters. Attributing this acceleration to a planetary mass companion, we 
demonstrate sensitivity to a planetary companion at 2 AU of mass 0.15 Mj. 

The residuals that we obtain from our fitting of 0.25 mas in each coordinate over 4.9 y indicate 
the feasibility of planet searching and characterization with radio astrometry of nearby low mass 
stars. These residuals can potentially be improved for this star and others through technical 
improvements to the data analysis including accounting for the motion of the star during the 
observing epoch and more exact removal of tropospheric fluctuations. The residuals indicate that 
astrometric jitter from stellar activity on scales of a stellar radius or larger is not a limiting factor 
in planet astrometry. Improvements in our analysis of systematic errors from stellar motion and 
tropospheric phase errors will permit us to explore sub-stellar-radius astrometric accuracy. 

Future RIPL papers will explore analysis of additional data from this star, improved analysis 
techniques, and analysis of other stars in the sample. Upgrades to the recording bandwidth of the 
VLBA (and GBT) will produce significant improvements in the capability to carry out observing 
programs of this kind. The improved sensitivity will allow more rapid detection of bright stars 
such as GJ 896A, thereby reducing any smearing effect from stellar motion, allow the use of closer 
calibrators for improved tropospheric corrections and reduced time on calibrator, provide higher 
signal to noise ratio detections of fainter stars, and enable us to detect a larger sample of stars. 

The National Radio Astronomy Observatory is a facility of the National Science Foundation 
operated under cooperative agreement by Associated Universities, Inc. This research has made 
use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion 



Laboratory, California Institute of Technology, under contract with the National Aeronautics and 
Space Administration. This research has made use of the SIMBAD database, operated at CDS, 
Strasbourg, France. This research was supported by NASA Origins grant NNX07AP19G. 
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Fig. 1. — Positions for GJ 896 A and primary and secondary calibrators. The distance from the 
primary calibrator, J2328+1929, to GJ 896A is 0.9 deg. The distance from the secondary calibrator, 
J2334+2010, to GJ 896A is 0.6 deg. The nearly linear arrangement and small separations make 
this an excellent set of calibrators for this star. 
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Fig. 2. — Initial model and astrometric data for GJ 896A. Top row shows astrometric data (black 
crosses) with the initial (optically-determined) astrometric model (green circles at the epoch of 
observation and red lines for the continuous calculation) . Positions in the top row are given in arcsec 
relative to a fiducial position of 23 ft 31 m , 19°56'. Bottom row shows residuals after subtraction of 
the model. The symbols in the top row are much larger than the errors; in the bottom row, errors 
include statistical and systematic contributions. 
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Final Model Astro 
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Fig. 3. — Fitted astrometric model and astrometric data for GJ 896A. Symbols are the same as in 
Figure El 



14 



Final Model Astro + Accel 
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Fig. 4. — Fitted astro-metric plus acceleration model and astrometric data for GJ 896A. Symbols 
are the same as in Figure [2 
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Fig. 5. — x 2 as a function of each parameter. The x-axis values are offsets from the best-fit value. 
The parameters explored are (1) right ascension (hours), (2) declination (degrees), (3) proper 
motion in right ascension (mas/y), (4) proper motion in declination (mas/y), (5) acceleration in 
right ascension (mas/y 2 ), (6) acceleration in declination (mas/y 2 ), (7) minimum time for the right 
ascension acceleration (days), (8) maximum time for the declination acceleration (days), and (9) 
parallax (mas). 
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Fig. 6. — Orbital motion for GJ 896B relative to GJ 896 A based on optical astrometric data. We 
plot the relative separation, velocity, and acceleration in the upper left, upper right, and lower 
left, respectiv ely. We plot in the lower right the historical data (black crosses) from which this 
orbit derived (JHeinta 1 1984 ) and the predictions of the orbital model (green circles at the time of 
measurement, red line for continuous time). We include an additional datum from 2007 that was 



not part of the fit (JMason et al. 



2001 



The plotted model curve extends from 1932 to 2017. 
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Fig. 7. — Periodogram power as a function of frequency for residuals after removal of astrometric 
and acceleration terms for the two coordinates. We also plot the 95% confidence threshold as the 
blue dashed line. The absence of a significant peak in the periodogram places limits on short-period 
companions. 
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Fig. 8. — Limits on companion mass as a function of semi-major axis for GJ 896A. Limits are derived 
from acceleration and angular separation terms. The space below the solid curve is allowed. 
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Table 1. Positions for Secondary Calibrator J2334+2010 



JD 



Flux 
(mJy) 



a 

(J2000) 



Aa 
(mas) 



5 

(J2000) 



AS 



mas 



2453820 


62.2 ±1.9 


23 34 14.1564960 


0.135 


20 10 28 


882640 


0.279 


2454482 


24.7 ±0.2 


23 34 14.1564897 


0.007 


20 10 28 


882867 


0.013 


2454585 


19.6 ±0.3 


23 34 14.1564998 


0.013 


20 10 28 


882864 


0.022 


2454680 


13.1 ±0.3 


23 34 14.1564977 


0.018 


20 10 28 


882824 


0.033 


2454793 


15.9 ±0.3 


23 34 14.1565001 


0.015 


20 10 28 


882793 


0.025 


2454978 


21.4 ±0.3 


23 34 14.1564972 


0.010 


20 10 28 


882838 


0.017 


2455101 


23.6 ±0.2 


23 34 14.1564926 


0.007 


20 10 28 


883053 


0.015 


2455226 


19.0 ±0.1 


23 34 14.1565001 


0.004 


20 10 28 


882733 


0.006 


2455467 


23.8 ±0.2 


23 34 14.1565155 


0.006 


20 10 28 


882853 


0.010 


2455591 


21.8 ±0.2 


23 34 14.1565109 


0.006 


20 10 28 


882888 


0.010 



Note. - - The first line is the mean position reported in Paper I. Errors in the position are errors 
from fitting and do not include systematic effects. 
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Table 2. Positions for GJ 896 A 



JD 



Flux 


a 


Aa 


S 


A5 


(A*Jy) 


(J2000) 


(mas) 


(J2000) 


(mas 





2453818.310 


929 ± 203 


23 31 52.4336850 


0.113 


19 56 13.712351 


0.171 


2453820.310 


1399 ± 298 


23 31 52.4342630 


0.183 


19 56 13.714736 


0.202 


2453821.310 


3895 ± 230 


23 31 52.4345400 


0.028 


19 56 13.715419 


0.045 



2454482.443 


328 ± 20 


23 31 52.4967356 


0.031 


19 56 13.570094 


0.081 


2454679.922 


278 ± 25 


23 31 52.5343729 


0.062 


19 56 13.704882 


0.128 


2454792.651 


2041 ± 52 


23 31 52.5301598 


0.015 


19 56 13.570517 


0.031 


2454978.109 


9512 ± 470 


23 31 52.5703944 


0.033 


19 56 13.606525 


0.060 


2455100.776 


1139 ±33 


23 31 52.5719845 


0.018 


19 56 13.599548 


0.051 


2455226.443 


138 ± 16 


23 31 52.5809398 


0.059 


19 56 13.443975 


0.101 


2455466.776 


2484 ± 88 


23 31 52.6126099 


0.021 


19 56 13.537846 


0.045 


2455591.443 


2018 ± 78 


23 31 52.6215506 


0.018 


19 56 13.383628 


0.044 



Note. - - The first three measurements are the results from Paper I. 



Table 3. Astrometric Parameters for GJ 896 A 



Parameter 


Unit 


Optical 


RIPL Astrometric 


RIPL Astrometric + Accel 


a 




23 31 52.17898 


23 31 52.179555 


23 31 52.179565 


Aa 


(mas) 


12.31 


0.409 


0.123 


5o 




19 56 14.1505 


19 56 14.141024 


19 56 14.139842 


A5 


(mas) 


9.51 


0.434 


0.136 


Ha 


(mas/y) 


554.64 ± 1.40 


571.333 ± 0.045 


571.169 ±0.014 


fJ-5 


(mas/y) 


-60.43 ± 1.08 


-60.781 ± 0.047 


-60.676 ± 0.014 


IT 


(mas) 


161.76 ±1.66 


159.864 ± 0.616 


160.010 ±0.182 


a a 


(mas/y 2 ) 






0.458 ± 0.032 


as 


(mas/y 2 ) 






0.087 ± 0.032 


4-a 


(JD) 






2454565 ± 32 


l 


(JD) 






2454879 ± 167 
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Table 4. Orbital Parameters for the GJ 896 AB System From Optical Astrometry 



Period 

(y) 


Semi-Major Axis 
(arcsec) 


Inclination 
(deg) 


n 

(deg) 


Periastron 

(y) 


Eccentricity 


CO 

(deg) 


359 


6.87 


123.5 


82.1 


2008 


0.20 


354.0 



